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1 Field Of The Invention 

2 The invention relates generally to the field of computer graphics, and more specifically to 

3 generation of smooth feature lines on subdivision surfaces representing surfaces of objects. 

4 Background Of The Invention 

^5 Generally, in computer graphics, objects are represented as surfaces, with the surfaces being 

0 

Q6 represented by meshes at a particular level. A mesh at a particular consists of a set of vertices, or 

01 

pj7 points in three-dimensional space, which are interconnected by edges. The edges define polygonal 

M 8 faces, which may be in the form of triangles, quadrilaterals, and so forth. In some computer graphic 

9 operations, it is desirable to generate a representation of a surface at a finer resolution than a current 

y 0 representation. There are several popular methodologies for generating a representation of a surface 

fU 1 at a finer resolution than a current representation, including a Catmull-Clark surface subdivision 

9 2 methodology, which is used in connection with surfaces defined by a quadrilateral mesh, and a Loop 

Q 3 surface subdivision methodology, which is used in connection with surface defined by a triangular 

1 4 mesh. Generally, both methodologies make use of respective subdivision rules at respective vertices 

1 5 defining the surface at a particular level in the mesh to generate a mesh in a next higher subdivision 

1 6 level. The surface of the respective object, which is referred to as the "subdivision surface" or "limit 

17 surface," is taken as being defined by a mesh as the subdivision level approaches infinity. 

1 8 A feature on a subdivision surface can be defined by a feature line in the mesh that defines 

19 the respective surface. A feature line can be in the form of a sharp crease or a smooth curve. For 

20 a smooth feature line, a normal vector, which is the vector that is perpendicular to a plane that is 

2 1 tangent to the surface, will vary continuously across the smooth feature line. On the other hand, for 

22 a sharp crease, the normal vector varies discontinuously across the crease, and in fact is not defined 

23 at the crease. However, a definition for a smooth curve can be derived from a definition for a sharp 

24 crease using one or more parameters that are used to define the sharpness of the curve across vertices 

25 in the mesh at a particular subdivision level. D. Zorin, "Stationary Subdivision And Multi- 

26 Resolution Surface Representation," Ph.D. Thesis, California Institute of Technology, Pasadena, C A, 
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1 1998, describes a methodology for generating a smooth feature line using parameters, but the 

2 methodology described there results in surfaces of relatively low quality, even in surfaces that are 

3 not smooth at some vertices in the surface topology. T. DeRose, et al., "Subdivision surfaces in 

4 character animation," SIGGRAPH 98 Conference Proceedings, Annual Conference Series, pages 

5 85-94, ACM SIGGRAPH, 1 998, describes a methodology for generating a smooth feature line based 

6 on applying a subdivision rule corresponding to sharp creases, as described in H. Hoppe, et al., 

7 "Piecewise Smooth Surface reconstruction," SIGGRAPH 94 Conference Proceedings, Annual 

.8 Conference Series, ACM SIGGRAPH, 1994, up to a selected fineness level, and a rule 

p* 

Q9 corresponding to smooth interior points thereafter. Applying the two distinctly different types of 

Q 

s$0 rules makes efficient evaluation of the resulting surface difficult. 

ru 

41 1 Summary Of The Invention 

iU 

2 

Jn2 The invention provides a new and improved system and method for generating a smooth 

Q3 feature line in a subdivision surface. 

4 In brief summary, the invention provides an arrangement for generating a representation of 

15 a feature in a surface defined by a mesh representation, the mesh comprising at a selected level a 

16 plurality of points including at least one point connected to a plurality of neighboring points by 

17 respective edges, the feature being defined in connection with the vertex and at least one of the 

18 neighboring points and the edge interconnecting the vertex and the at least one of the neighboring 

19 points in the mesh representation. The feature generating arrangement comprises a weight vector 

20 generator module and a feature representation generator module. The weight vector generator 

21 module is configured to generate at least one weight vector based on a parameterized subdivision 

22 rule defined at a plurality of levels, for which a value of at least one parameter differs at at least two 

23 levels in the mesh. The feature representation generator module configured to use the at least one 

24 weight vector and positions of the vertex and the neighboring points to generate the representation 

25 of the feature. 

26 One embodiment of the invention comprises an arrangement for generating a smooth feature 

27 line in the subdivision surface. 

i 
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1 Brief Description Of The Drawings 

2 This invention is pointed out with particularity in the appended claims. The above and 

3 further advantages of this invention may be better understood by referring to the following 

4 description taken in conjunction with the accompanying drawings, in which: 

5 FIG. 1 depicts a computer graphics system including an arrangement for generating a smooth 

6 feature line in a subdivision surface, constructed in accordance with the invention; 

M 

^7 FIG. 2 depicts a mesh representing a surface to which the Loop surface subdivision 

01 8 methodology has been applied; 

fU 

9 FIG. 3 depicts a stencil useful in understanding the Loop surface subdivision methodology; 

ru 

» 10 FIGS. 4 A through 4E depict a mesh representing a surface to which the Catmull-Clark 

|j 1 surface subdivision methodology has been applied; 

o 

H2 FIG. 5 depicts a stencil useful in understanding the Catmull-Clark surface subdivision 

1^3 methodology; 

1 4 FIG. 6 depicts a stencil useful in understanding a methodology for generating a sharp crease 

15 on a subdivision surface; 

16 FIG. 7 depicts a stencil useful in understanding a methodology for generating a smooth 

1 7 feature line used in connection with subdivision surfaces defined by triangular meshes; 

1 8 FIGS. 8 and 9 depict tables of coefficients used in the methodology used in connection with 

19 surfaces defined by triangular meshes described in connection with FIG. 7; 

20 FIG. 10 depicts a flow chart describing operations performed by the computer graphics 

2 1 system in connection with the methodology used in connection with surfaces defined by triangular 

22 meshes described in connection with FIG. 7; 

23 FIG. 1 1 depicts a stencil useful in understanding a methodology for generating a smooth 

24 feature line used in connection with subdivision surfaces defined by quadrilateral meshes; 
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1 FIGS. 12 and 13 depict tables of coefficients used in the methodology used in connection 

2 with surfaces defined by quadrilateral meshes described in connection with FIG. 1 1 . 

3 Detailed Description of an Illustrative Embodiment 

4 FIG. 1 depicts a computer graphics system 10 including an arrangement for generating a 

5 smooth feature line in connection with a subdivision surface, constructed in accordance with the 

6 invention. With reference to FIG. 1 , the computer graphics system includes a processor module 1 1 , 
g7 one or more operator input devices 12 and one or more display devices 13. The display device(s) 
^8 13 will typically comprise a frame buffer, video display terminal or the like, which will display 
fij9 information in textual and/or graphical form on a display screen to the operator. The operator input 
*J0 devices 1 2 for a computer graphics system 1 0 will typically include a pen 1 4 which is typically used 
^1 in conjunction with a digitizing tablet 15, and a trackball or mouse device 16. Generally, the pen 
f|2 14 and digitizing tablet will be used by the operator in several modes. In one mode, the pen 14 and 
L^3 digitizing tablet are used to provide updated shading information to the computer graphics system. 
H 4 In other modes, the pen and digitizing tablet are used by the operator to input conventional computer 

a 

5 graphics information, such as line drawing for, for example, surface trimming and other information, 

16 to the computer graphics system 10, thereby to enable the system 10 to perform conventional 

17 computer graphics operations. The trackball or mouse device 16 can be used to move a cursor or 

1 8 pointer over the screen to particular points in the image at which the operator can provide input with 

19 the pen and digitizing tablet. The computer graphics system 10 may also include a keyboard (not 

20 shown) which the operator can use to provide textual input to the system 10. 

21 The processor module 1 1 generally includes a processor, which may be in the form of one 

22 or more microprocessors, a main memory, and will generally include one a mass storage subsystem 

23 including one or more disk storage devices. The memory and disk storage devices will generally 

24 store data and programs (collectively, "information") to be processed by the processor, and will store 

25 processed data which has been generated by the processor. The processor module includes 

26 connections to the operator input device(s) 12 and the display device(s) 13, and will receive 

27 information input by the operator through the operator input device(s) 12, process the input 

28 information, store the processed information in the memory and/or mass storage subsystem. In 
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1 addition, the processor module can provide video display information, which can form part of the 

2 information obtained from the memory and disk storage device as well as processed data generated 

3 thereby, to the display device(s) for display to the operator. The processor module 1 1 may also 

4 include connections (not shown) to hardcopy output devices such as printers for facilitating the 

5 generation of hardcopy output, modems and/or network interfaces (also not shown) for connecting 

6 the system 10 to the public telephony system and/or in a computer network for facilitating the 

7 transfer of information, and the like. 

1^8 The invention provides an arrangement for generating a smooth feature line in a three- 

Q9 dimensional subdivision surface. The subdivision surface, in turn, is defined by control points that 

U 

OHO are interconnected by edges to form a three-dimensional mesh. Generally, a subdivision surface is 

jjj 1 initially defined by a mesh at a particular degree of granularity or fineness. Using one of several 

^2 methodologies, the mesh can be refined through a series of levels of increasing subdivision levels, 

9 13 with the subdivision surface being the limit as the subdivision level approaches infinity. In the 

4 following, the each mesh level will be identified an index, with one mesh level being identified by 

□ 5 index "j" and the next higher subdivision level mesh being identified by index "j+1." The 

^ 6 subdivision surface is essentially defined by the positions of the surface's "limit points," which are 
the essentially the control points as the subdivision level approaches infinity, and the orientations 

1 8 of the associated normal vectors. Each normal vector defines the orientation of a tangent plane at 

19 the respective point on the subdivision surface and, in turn, can be defined by the cross product of 

20 vectors that are tangent to the subdivision surface at the respective point. 

21 Before describing the operations performed by the smooth feature line generating 

22 arrangement, it would be helpful to describe two surface subdivision methodologies, namely the 

23 aforementioned Loop surface subdivision methodology and a Catmull-Clark surface subdivision 

24 methodology, with which the arrangement is used.. 

25 Loop's surface subdivision methodology will be described in connection with FIGS. 2 and 

26 3. FIG. 2 depicts a portion of an illustrative surface defined by a mesh, and a subdivision surface 

27 derived therefrom, and FIG. 3 depicts a stencil useful in understanding the surface subdivision 

28 methodology. Generally, in Loop's surface subdivision methodology, each triangular face in the 

29 original mesh is split into a plurality of subfaces, the subfaces defining the next higher subdivision 
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level mesh. The vertices of the next higher subdivision level mesh are positioned using weighted 
averages of the positions of the vertices in the original mesh. More specifically, and with reference 
to FIG. 2, a mesh 20 at a "j" subdivision level is depicted that includes a vertex v q (vertex v q =v q (0)) 
21(0) located at position c^q), and a plurality "K" of surrounding points v q (l) through v q (K) 
(generally, v q (k)) 21(1) through 21(K) (generally identified by reference numeral 21(k)) located at 
positions at respective positions c*(l) through d(K), the points v q (k) 21(k) being those points in the 
mesh 20 that are connected to vertex v q by edges. ("K," the number of points connected to vertex 
v q will sometimes be referred to as the vertex's "valence.") The stencil depicted in FIG. 3 is useful 
in understanding the ordering of the index "k" for the points surrounding the respective vertex v q . 
Denoting the set of the indices of the points v q (k) in the mesh that are connected to vertex v q by 
N(qj), a mesh at the next higher subdivision level "j+l" corresponds to vertex v* q 22(0) located at 
position d +1 (q) and a set of surrounding points v q (l) 22(1) at located at positions and connected 
to the vertex v q 22(0). The higher subdivision level mesh is constructed by providing the vertex v' q 
22(0) at a position d +l (q) that is determined by 




(1) 



and points v' q (l) 22(1) at positions c* +1 (l) that are determined by 



* y+, (/) = l[c J (q)+c J (k)] + ^[c J (k-l)+c J (k+l)],k=\,...,K (2), 



where the weighting factor a(K) is given by 




(3). 
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1 Equations (1) and (2) can be summarized in a single equation defining a subdivision rule using a 

2 subdivision matrix S T K (where index "T" refers to "Triangular") as follows 

3 c J+l = S TtK c j (4), 



which may be written component-wise as 



(5), 



nj6 



ni 
a 

b 
ru 



where 



\-a(K) ifl=0,m=0 
a(K) 

— — if I = 0,m= \,...,K 



— { 



3 
8 

J_ 
8 



0 



«//= 1,. ,.,K,m= 0 
or I = m = \,...,K 

ifl=2,...,K,m=l-\ 
or I = \,m = K 
or 1= l,...,K, m = /+ 1 
or/= K,m= 1 



(6), 



8 where indices "1" and "m" range from "zero" to, and including, "K." 

9 Equations ( 1 ) and (2), or equivalently equations (4) through (6), are applied by selecting each 

10 point in the mesh at subdivision level "j" as a vertex, to provide points for the mesh at the 

11 subdivision level "j+1." Loop's surface subdivision methodology can be applied recursively to 

1 2 provide a mesh at any desired subdivision level. It will be appreciated that equation (2) will provide 
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1 that the point v' q (l) 22(1) that is created when the methodology is applied to point v q 21(0) as the 

2 vertex will be in the same position when the methodology is applied to point v q (k) 21(k) as the 

3 vertex. 

4 It will be appreciated that, the Loop surface subdivision methodology provides 

5 (i) one new point (reference points v f q (l) described above), at location c j+1 (1) in the next higher 

6 ("j+l") subdivision level mesh, for, and generally somewhere between, each pair of points in the 

7 level f, j" mesh, and 

08 (ii) a relocated point (reference vertex v' q described above), at location c* +1 (q) in the "j+l" 

019 level mesh for, and generally somewhat near, each point in the "j" level mesh, 

ru 

'!l0 which are interconnected to form the "j+l" level mesh. 

m 

a 1 1 The Catmull-Clark subdivision methodology produces smooth surfaces using a small number 

2 of neighboring vertices. The Catmull-Clark surface subdivision methodology will be described in 

&3 connection with FIGS. 4 A through 4E and 5. FIG. 4 A depicts an illustrative surface defined by a 

jolt 

U4 mesh, and FIGS. 4B through 4E are useful in understanding operations performed in generating a 

1 4 5 subdivision surface derived from the surface depicted in FIG. 4A. FIG. 5 depicts a stencil useful in 

1 6 understanding the surface subdivision methodology, and in particular is useful in understanding the 

1 7 ordering of the index "k" for the points surrounding the vertex v q (vertex v q =v q (0)) around which the 

1 8 surface subdivision methodology is to be performed. With reference to FIG. 4A, that FIG. depicts 

19 amesh 30 consisting of four quadrilaterals 31(1) through 31(4). Each quadrilateral is referred to as 

20 a face. The quadrilaterals are defined by points 32(0), which is common to all of the quadrilaterals 

21 31(1) through 31(4), and other points 32(1) through 32(8). Generally, the Catmull-Clark surface 

22 subdivision methodology is performed in a series of iterations, including 

23 (i) a face point generation iteration, in which, for each face, a face point 33(1) through 33(4) 

24 is generated as the average of the points 32(p) defining the respective face (reference FIG. 4B); 

25 (ii) an edge point generation iteration in which, for each edge, a new edge point 34(1) 

26 through 34(4) is generated as the average of the midpoint of the original edge with the two new face 
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1 points of the faces adjacent to the edge (reference FIG. 4C; the face points generated during the face 

2 point generation iteration are shown as circles); 

3 (iii) a vertex point iteration in which a vertex point 3 5 is generated in relation to the positions 

4 of the new face points generated in iteration (i) for faces that are adj acent to the original vertex point, 

5 the positions of the midpoints of the original edges that are incident on the original vertex point, the 

6 position of the original vertex point, and the valence of the original vertex point (reference FIG. 4D; 

7 the face and edge points generated during the face point and edge point generation iterations are 
1^8 shown as circles); and 

a 

Q9 (iv) a mesh connection step (reference FIG. 4E) in which: 

01 

jjj 0 (a) each new face point is connected to the edge points of the edges defining the original 

*H 1 face; and 

fU 

%]2 (b) each new vertex point is connected to the new edge points of all original edges 

^ 3 incident on the original vertex point. 

f44 More specifically, in the Catmull-Clark subdivision methodology, denoting, for a vertex v q 
at location d(q), 

16 (i) the set of indices of first-order neighboring points that are connected thereto (for 

17 example, for point 32(0), point 32(2) and 32(8)) by N e (q,j), and 

1 8 (ii) the set of indices of second-order neighboring points that are opposite to vertex d(q) with 

19 respect to the level "j" faces (for example, for point 32(0), point 32(1)) that are incident with vertex 

20 c^byN^qj) 

21 a mesh at the next higher subdivision level "j+l" is constructed as follows. In the face point 

22 generation iteration, the face points are generated and located at positions c i+1 (m i ) determined as 

23 follows 
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= \ (c>{q) + <f(/,) + *>(/ l+I ) + c>(4 



m i € N f (q,j+ l),/„/ l+I € N e (qj),n e N f (qj) 



2 In the edge point generation iteration, the edge points are generated and located at positions c^'Oj) 

3 determined as follows: 



4 e N e {q,j+ l),*, e N e (q,j\m i _ x ,m i e N f (q,j+ l) 



j 5 In the vertex point generation iteration, the new vertex points are generated and located at positions 
' 6 as follows: 



8 In terms of only points d, that is substituting for the last term in equation (9), 



(10). 
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1 Generally, for quadrilateral meshes, the arrangement makes use of the Catmull-Clark 

2 methodology, except for the case of K=3 in equation (9). In that case, the arrangement makes use 

3 of 

c J+x {q) = (l-8rM*)+4r-jp Z^(/)+4 r ^ 2> +, (m) (11), 
3 

5 where y = — . Equation (1 1 ) would correspond to the Catmull-Clark methodology (equation (9)) 

{=* 3 o 

a 
a 

81 1 3 1 

pj6 with y - — — , except that, for K=3, / = — instead of — . In terms of only points c 1 (compare 
! 4K 38 12 

H 

€1 

fU7 equation (10)) 



a 

w 
a 

M 



c J+l (q) = (l-7y)c J (g)+6y^ !>(/)+/ ^ J]c J (m) ( i 2 ). 



nJ 9 These equations can be summarized by a single equation defining a subdivision rule using the 

10 subdivision matrix S Q K (where index "Q" refers to "Quadrilateral") as follows: 

11 c J+} = S Q , K c J (13), 



1 2 which may be written component- wise as 
13 



14 where 
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1 4K 

1JL 
1 



* 2 

3 
8 

J_ 
16 



— < 



4 
0 



i/ / = 0,m = 0 

ifl= 0,m= l 9 ... 9 K 

if 1 = 0 9 m = K+ 1,... 9 2K 

if /= l 9 ... 9 K,m = Oor I 

ifl= l,m= K 

or /= 2,...,K,m= I- 1 

or / = l,...,i£- 1, m = /+ 1 

or / = K 9 m = 1 

or I - 1, ra = 2 AT 

or/=2,...,,K, ra=AT+/-l 

or / = l 9 ... 9 K 9 m = 1 

ifl= K+l 9 ... 9 2K 9 m = 0,/- KJ- K+l,orl 

orl=2K 9 m=0 9 K 9 \or 2K 

otherwise 



(15), 



where indices "1" and "m" range from "zero" to, and including, "2K." 

As is apparent from the above discussion, the Catmull-Clark surface subdivision 
methodology applied to a quadrilateral mesh at one "j-th" subdivision level will produce a 
quadrilateral mesh at the next higher "j+l-st" subdivision level having the same number of points 
as the mesh at the "j-th" subdivision level, but the points in the higher ("j+1 ") subdivision level may 
be at different positions than in the "j-th" level. Equations (7) through (12), or equivalently 
equations (13) through (15), are applied by selecting each point in the mesh at subdivision level "j" 
as a vertex, to provide points for the mesh at subdivision level The Catmull-Clark surface 
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subdivision methodology can be applied recursively to provide a mesh at any desired subdivision 
level. 



Sharp creases in a subdivision surface may be generated by local modifications of the surface 
subdivision equations described above (equations (6) through (8) in the case of a triangular mesh, 
equations (13) through (15) in the case of a quadrilateral mesh). This will be illustrated in 
connection with FIG. 6, which depicts a stencil that is useful in understanding the ordering of the 
index "k" for the points on opposing sides of a respective vertex v q (vertex v q =v q (0)) through which 
a crease is to occur. By marking the edges of the mesh at which a crease is to occur, chains of edges 
can be defined that will define the position of a crease on a subdivision surface. As in the cases of 
both triangular and quadrilateral meshes, the positions d +1 of vertices and points in the higher 
"j+l-st" subdivision level can be determined from positions of vertices and points d in the "j-th" 
level using a defining a subdivision rule comprising a subdivision matrix S c (where index "C" refers 
to "crease") as follows: 




(16), 



which may be written component-wise as 




(17), 



m 



where 



3 

4 
1 



ifl=0,m=0 





(18), 



2 
0 



ifl=lor2,m=0orl 
otherwise 
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1 where indices "1" and "m" range from "zero" to, and including, "two." This subdivision rule can be 

2 used for creases in the cases subdivision surfaces defined by both triangular meshes and quadrilateral 

3 meshes. 

4 The invention provides a system and method for generating smooth feature lines in 

5 subdivision surfaces. As with a sharp crease, a smooth feature line is defined by marking a series 

6 of edges at a given subdivision level, which edges are to define the feature line. The lowest 

7 subdivision level at which the edges that are associated with a smooth feature line are marked will 

M 8 be referred to as the definition level "j D " for the smooth feature line In addition, a sharpness 

Q 

q 9 parameter will be associated with each marked edge in the definition level j D for the smooth feature 

j^j 0 line. In one embodiment, the value of the sharpness parameter is in the form of a real number in the 

K 1 interval [0, 1 ], that is, in the interval between "zero" and "one," inclusive of the endpoints. If an edge 
has a sharpness parameter with the value "zero," there is no crease along the edge. On the other 

* 13 hand, if an edge has a sharpness parameter with the value "one," the edge forms part of a sharp 

PJ14 crease as described above. 

a 

jj5 Generally, a smooth feature is obtained by using a subdivision rule that generalizes the 

fyi6 subdivision rules described above (reference equations (4) through (6) in the case of a triangular 

17 mesh and equations (13) through (15) in the case of a quadrilateral mesh). The subdivision rule 

18 makes use of two parameters s, and s 2 which are the respective values of the sharpness parameter 

19 on two sides of a vertex across which a smooth feature line is to be generated. A subdivision 

20 surfaced limit surface, that is, the surface defined by the subdivision surface's limit points, is smooth 

21 in the neighborhood of a smooth feature line. That is, for each point on the limit surface that is 

22 associated with a smooth feature line, that is, for each point for which the values of the sharpness 

23 parameters on respective sides of the point are other than "one," the orientation of the vector that is 

24 normal to the limit surface will vary continuously. However, for the limit surface at a point for 

25 which the sharpness parameter have the value is "one," which will be the case at a sharp crease, the 

26 normal vector will not be defined, since the orientation of the normal vector on respective sides of 

27 the point will differ. 

28 Operations performed by the computer graphics system 1 0 in connection with generating a 

29 smooth feature line will be described in connection with FIGS. 7 through 10, in the case of a 
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subdivision surface defined by a triangular mesh, and in connection with FIGS. 1 1 through 14, in 
the case of a subdivision surface defined by a quadrilateral mesh. In both cases, the smooth feature 
line will be defined by the positions of limit points, as well as the normal vectors at the respective 
limit point, in a region of a subdivision surface that is proximate the line corresponding to the edges 
that have been marked. As noted above, the normal vectors correspond to the cross product of two 
tangent vectors at the respective limit point, namely, a tangent vector along the smooth feature line 
at the respective limit point and a tangent vector across the smooth feature line at the respective limit 
point. 

With reference to FIG. 7, that FIG. depicts a portion of a subdivision surface defined by a 
triangular mesh with vertex v q (vertex v q =v q (0)) and points v q (l) through v q (K) with indices 
numbered as depicted in the FIG. The vertices are labeled such that a crease, which is to define part 
of a smooth feature line, passes through point v q ( 1 ), vertex v q and point v q (L+ 1 ) along edges v q ( 1 ), v q 
(which will be referred to as edge (1,0)) and v q ,v q (L+l) (which will be referred to as edge (0,L+1)). 
The value of the sharpness parameter associated with the edge (1,0) will be referred to as and the 
value of the sharpness parameter associated with the edge (0,L+1) will be referred to as s 2 . Since 
edges (1,0) and (0,L+1) have been indicated as being edges that define a crease, the value of the 
sharpness parameter will be non-zero. The values of the sharpness parameters associated with the 
other edges that are on the vertex v q may be"zero," or non-zero. Generally the subdivision rule 
described in equation (4) (or, equivalently, equation (5)) can be used to determine the positions of 
respective vertices, using the subdivision matrix 
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(l- s 3 ){l-a(K))+-s 2 ifl = 0,m= 0 
/ \a(K) 1 

(l-5 3 j-^-+-5 3 if 1 = 0,m= lorL+1 

(l-^s)^ tfl=0,m=2,...,L 

orl=0,m = L+2,...,K 

3 1 

-+-5 2 ifl=\,m=Oor\ 
3 1 

+ if I = L+ \,m= Oor L+ I 



8 8 1 

— (l-j 2 ) if I = \,m = lor K 
8 

^(l-i,) if 1= L*\,m= Lor L + 2 

3 (19) ' 
8 if l=2,...,L,m= 0 

or 1 = L+ 2,...,K,m= 0 
or l- m - 2,...,L 
or I - m - L+2,...,K 

^ if 1= 2,...,L,m= l-l 

or 1= L+ 2,...,K, m- I- 1 
or /= 2,...,L, m= 1+ 1 
orl= L+ 2,...,K,m= 1+ 1 
or / = A', w = 1 

0 otherwise 
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1 where, as in equation (6), indices "1" and "m M range from "zero" to, and including, "K," and where 

2 s 3 = — (jj + s 2 ) , the average of the values of the parameters s, and s 2 associated with the edges 

2 

3 (l,0)and(0,L+l). 

4 To ensure that the directions of the tangent vectors vary continuously at and near irregular 

5 vertices, the value(s) of the sharpness parameter(s) are gradually reduced as the level of the mesh 

6 increases. In the case of smooth feature lines, the repeated update of the value of the sharpness 

u 

q7 parameter as the level of the mesh increases would provide a value that approaches zero, which, in 

^8 turn, will provide that the directions of the tangent vectors vary continuously at the limit surface. 

fU 9 In one embodiment, for a vertex v q for which the sharpness parameters s, and s 2 have the same value 

yjO "s" on both sides of the vertex (that is, s,=s 2 ), the sharpness parameter update function is 

ITJ 

all 

HI 



18 



s{ J+ l) = (s{j)f (20), 



CJ2 that is, s(j+l), the value of the sharpness parameters s, and s 2 at the higher M j+l-st" subdivision level 

M 

' 1 3 mesh, is the square of s(j), the value of the sharpness parameters s x and s 2 at the lower ("j-the") level 

14 mesh. On the other hand, for a vertex v q for which values of sharpness parameters Sj and s 2 are not 

15 constant across the vertex (that is, s, *s 2 ), the sharpness parameter update functions are a function 

16 of the values of both sharpness parameters s, and s 2 . In one embodiment, the function is selected to 

17 be 



.(/+») - (f*,o)^o)) 2 (2l) - 



19 and 
20 



2 1 essentially providing a linear combination corresponding to quadratic splines. 
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Since the values of the sharpness parameters Sj and s 2 vary from level to level, the smooth 
curve subdivision matrix (S sc T KL (s„s 2 )) Im (equation (19)) will also not be constant from level to 
level. Accordingly eigen-analysis, which is conventionally used to determine positions of limit 
points and orientations of tangent vectors for subdivision surfaces, is not used here. As suggested 
above, a subdivision surface is defined by a collection limit points and normal vectors (or, 
equivalently, two tangent vectors for each normal vector), which, in turn, are provided by iteratively 
applying equation (4) from a definition level M j D " for the mesh to the level j=«>. The definition level 
j D is the lowest subdivision level "j" for which the respective smooth feature line is defined in the 
subdivision surface. Accordingly, in order to determine the positions of the limit points, the infinite 
subdivision matrix 



is evaluated, where S sc T K L (s 1 (j),s 2 (j)) is the subdivision matrix described above in connection with 
equation (19) for sharpness parameters corresponding to the "j-th" level in the matrix product on the 
right-hand side of equation (23). For S scTKLL p(s„s 2 ) on the left-hand side of equation (23), 
arguments s x and s 2 refer to the sharpness parameters at the definition level of the smooth feature 
line, that is, s^O'd) and s 2 =s 2 (j D ), and the subscript "LP" refers to "Limit Point." Since matrix 
multiplication is not commutative, the order of the factors in the matrix product in equation (23) is 
important. In equation (23), the order of the factors from left to right in the product extends from 
j=oo on the left to j=j D on the right, as suggested by equation (4). 

The matrix product (equation (23)), converges to a "K+l"-by-"K+r matrix that has "K+l" 
rows, with the components comprising each row being identical to the components of the other rows, 
in the same order. A vector of weight values 1 LP for determining the position a(q) of the limit point 
is obtained as follows: 





(24), 
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where vector v LP is a vector of weights that are used in determining the position of a limit point of 
a subdivision surface using the Loop subdivision methodology in the absence of a smooth feature 
line, that is, 



6){K) 1 1 1 



ni(K)+ K' m(K)+ K' m(K)+ K' 'm(K)+ K) 



(25), 



where 



3K 

®( K ) = TTT^ (26), 



Sa(K) 



where a(K) is as defined in equation (3). The position a(q) of a limit point associated with a vertex 
v q (0) is computed from the position c jD (q) of the vertex v q (0) and the positions c jD (k) ,k=l,...,K, 

of the neighboring points v q (k) around the vertex v q (0), in the mesh corresponding to the definition 
level j D of the smooth feature line in the subdivision surface, or to any subdivision level higher than 
the definition level, using the components of 1 LP as weight values. The components of limit point 
weight vector l LP (s,,s 2 ) can be satisfactorily approximated by second degree polynomials in s x and 
s 2 using six values per component and per valence. 

K 

If L = — , symmetry considerations allow the number of values per component to be 

reduced. In that case, the components of limit point weight vector 1 LP satisfy the following symmetry 
relations: 

1. The "i-th" component of limit point weight vector 1 LP equals the "K+2-i-th" 

K 

component of 1 LP , that is, [l LP ) = [l Lp ) > for i = — + 2, . . . , K . This is a 
result of the fact that the subdivision matrix (equation (19) is invariant upon 
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reflection with respect to the crease line. As a result, approximations need only be 

K 

generated for components / = 0,1, . . . , — + 1 . 

2 



Further, since the subdivision matrix (equation (19)) is also symmetric in s, and s 2 , 
(^C 5 !'^)). ~ (^(^'^i))^ 2 > and so approximations need only be generated 



for i = 0,1,... , 



K 








, where 11 









" refers to the greatest integer in the quotient 



x 
it .it 

y 



Further, since, if the particular values of the parameters Sj and s 2 were interchanged, 
the vertex v q would behave in the same way, {l Lp {s x , s 2 )) q = [l LP {s 2 , s x )^ , and 

so the approximation may be constructed using only symmetric, second-degree 
polynomials in Sj and s 2 , that is, using four coefficients; and 

Further, if "K" is a multiple of "four," then, from points (1) and (2) above, 

{^Lp{ s \9 s 2f) K = {Jlp{ s 2> s i})k • Accordingly, these components can be 

approximated using symmetric second-degree polynomials using four coefficients. 
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K 

In one embodiment, regardless of whether L = — , the approximating polynomials are computed 

using a least-squares Chebyshev approximation methodology. The matrix product (equation (23)) 
is evaluated at the points of a grid in the (s„s 2 ) region made up of roots of Chebyshev polynomials 



{s x ,s 2 ) = 



cos 







{{■ '1 






{ J+ 2) 




,cos 


N 


N 





71 



(27), 



forij=0,...N-l. In one embodiment, "N" is taken to be "eight," so that there are "sixty- four" sample 
points in equation (27). The coefficients \ (j=0,...,3 in the symmetric case, or j=0,...,5 in the 
asymmetric case) for the polynomial 



(28), 



K 

in the symmetric case ( / = Oor — + 1 ), or the polynomial 

4 



( 1 lp\ * 6,0 + ^ + ^2 + b n s 2 x +b iA sl + b i5 s x s 2 



(29), 



in the asymmetric case ("i" otherwise), are determined by a least squares methodology. In that 
operation, the computer graphics system 10 computes the matrix product (equation (23) evaluated 
at points (s,,s 2 ) as defined in equation (27)). Thereafter, a set of N 2 polynomials (equation (28) 
and/or (29)) are generated using the values of Sj and s 2 at the respective points (equation (27)) in the 
(s l5 s 2 ) region, with each polynomial being equal to the respective "i-th" component of the row of the 
matrix product. The computer graphics system 10 then determines the values of the coefficients by 
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in the polynomials using the least squares approximation methodology. After the values of the 
coefficients have been determined, they are used to generate the values of the components (l LP ) i 
of limit point weight vector 1 LP using equations (28) and (29). 

The error in the least squares approximation has an oscillatory behavior over the region (s { ,s 2 ) 
and will be minimized if the amplitude of the oscillations is distributed evenly over the region. This 
can be accomplished by relaxing the constraints imposed by a weighted least squares methodology 
for points that are well-approximated and tightening the constraints for the points that are poorly 
approximated. In the weighted least squares methodology, after the values of the coefficients b^ 
have been evaluated using the least squares methodology described above, the polynomial (equation 
(28) and/or (29)) is evaluated at each point in the (s l9 s 2 ) region, the value compared to the value 
generated for the respective component of the matrix product (equation (23)), and a weight value 
determined that reflects the difference therebetween. Thereafter, the weight values can be used in 
constraining the approximations of the values of the coefficients b {j in an uneven manner. 

FIG. 8 depicts a table of coefficients by for K=4, 6,..., 16, with the indexing, in the form of 



the C programming language, b 



K 

y- 2 



[/][ y] , where "i" is the index of the component of limit 



point weight vector 1 LP , taking values in the range / = 0, . . . , 



K 

4 + 1 



, and "j" is the index of the 



coefficient b {j of the approximating polynomial. 

After the limit point weight vector 1 LP associated with a vertex v q (0) has been generated as 
described above, it can be used along with the position d(q) of the vertex v q (0) and positions d(k) 
of the points v q (k), k=l,...,K neighboring the vertex v q (0) to determine the position a(q) of the limit 
point associated with the vertex v q (0) as follows 



= Z { l Lp{^ 2 ))c J (i) 

i=0 



(30), 
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where, as suggested above, i e qu Nyq,j) (reference FIG. 7), and where (l^s,^))! refers to the 



"i-th" component of l LP (s„s 2 ), which, in turn, corresponds to the limit point weight vector 1 LP for 
sharpness parameter values s { and s 2 at the definition level j=j D for the smooth feature line or at any 
subdivision level higher than the definition level j D , and indices 0,...,K are as shown in FIG. 7. 

One embodiment makes use of a similar polynomial approximation methodology to obtain 
good approximations for vectors \q and l s of weight values that are to be used in generating the 
tangent vectors, tangent vector weight vector l c being used in generating the tangent vector along 
the smooth feature line and tangent vector weight vector l s being used in generating the tangent 
vector across the smooth feature line at the limit point associated with the vertex v q (0) with which 
the respective tangent vector weight vectors is associated. In that methodology, a vector 1 C (J) is 
generated as follows: 



where the two non-zero components of the row vector on the right hand side are a "one" at position 

K 

"one" in the row vector, and a "negative one" at position — + 1 , and "•" refers to the dot product. 

After generating vector 1 C (J) according to equation (31), the value is normalized and the limit taken 

as J-°o; 



1 C (J) = (0X0 9 ...rW,..)>f[s Kj j J .(s l Q) t 3 2 (j)) 



(31), 



l c = lim 



1 C {J) 



(32), 



1 C (J)\ 
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1 where ||v|| = ^f^v y 2 , that is, the Euclidean norm. The components of vector l c (equation (32)) 

2 form a sequence of weight values that are used in multiplying the respective positions of the vertex 

3 v q (0) and the points v q (k), k=l ,...,K, around the vertex v q (0), as shown in FIG. 7, the sum of which 

4 defines the tangent vector along the smooth feature line at the point on the subdivision surface 

5 associated with the vertex v q (0). 

6 As with the components of limit point weight vector 1 LP , the components of tangent vector 

M 

0? weight vector l c can be approximated by polynomials in s x and s 2 , in this case third degree 

§8 polynomials. Accordingly, the approximations of the components of tangent vector weight vector 

fU9 l c can be efficiently generated by using up to ten numbers per component and per valence. If 

M 

pjO L - — , symmetry considerations allow the number of polynomial coefficients that are used in 

a 

FUl generating components of tangent vector weight vector l c to be reduced. In that case, the 

Q 2 components of tangent vector weight vector l c satisfy the following symmetry relations: 

Q 

Hi 3 1 . The "i-th" component and the M k+2-i-th" component of tangent vector weight vector 

14 l c are equal, that is, (l c ) = [l c )^ +2 . , for i = — + 2, . . . , K . This is a result of the 

1 5 invariance of the subdivision matrix (equation (19)) upon reflection with respect to 

1 6 the crease line. As a result, approximations need only be generated for components 

K 

17 i = 0,l,...,y+l; 



18 2. Further, since (/ c (^,5 2 )) = - (l c ( s 2> s i)) K +2 .» approximations need only be 



19 generated for components / = 0,1, • ■ • , 
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1 3 . Further, since, if the particular values of the parameters s, and s 2 were interchanged, 

2 the vertex v q would behave in the same manner except that the orientation of the 

3 tangent vector would be reversed, (j c [s x , s 2 = - [l c (s 2 , s x )) q , in which case 

4 the approximation may be constructed using only anti-symmetric, third-degree 

5 polynomials in s, and s 2 , using four coefficients; and 

6 4. Further, if K is a multiple of "four," then, from points (1) and (2) immediately above, 

M 

07 (jc( 5 i> 5 2))* = " (^c( 5 2> 5 i))* • Accordingly, these components can be 

flfl 
ry 

approximated using anti-symmetric polynomials using four coefficients. 

I1J9 In one embodiment, the approximating polynomials are also computed using a least squares 

2 

f|0 Chebyshev approximation methodology similar to that used in connection with the limit point 

jjj 1 approximation described above. The matrix product of equation (3 1) is evaluated at points of a grid 

M2 in the (s„s 2 ) region made up of roots of Chebyshev polynomials (equation (27)). The coefficients 

ejp bjj (j=0,...,3, in the anti-symmetric case, or j=0,...,9, in the non-symmetric case) for the polynomial 



14 



16 



K 

15 in the anti-symmetric case ( / = 0 or — + 1 ), or the polynomial 

4 



(0, = b io + b n s i + b i2 s 2 + &, 3 J? + ft| 4 Ji*2 + 

b i5 S 2 + */6 J ? + b il S \ S 2 + */8 J l ff 2 + */9 J 2 



(34), 



17 in the non-symmetric case ("i" otherwise), are determined using the least squares methodology as 

18 described above. 
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In a manner similar to the generation of the tangent vector along the smooth feature line 
(reference equations (31) through (34) above), a vector l s of weight values is used in generating the 
tangent vector across the smooth feature line. The components of tangent vector weight vector l s 
comprise a sequence of weight values that are used to multiply the positions of the vertex v q (0) and 
points v q (l) through v q (K) around vertex v q (0) (reference FIG. 7), the sum of which comprises the 
tangent vector across the smooth feature line. In that methodology, a vector 1 S (J) is generated as 
follows: 

f . 2*(0) . 2/r(l) . 2*(K-l)) A / M M x 

l s yJ) = [0,sin— —,sm— — ,...,sin J [ [ S^r^s^jj^j)) (35). 

Vector 1 S (J) is normalized and the limit taken as J-~, in the same manner as described above in 
connection with equation (32), to provide tangent vector weight vector l s . 

K 

When L = — , tangent vector weight vector l s is 



f . 2*(0) . 2n(\) . 24*- 01 
l s = 10, sin — ,sin— — ,...,sin — J (36), 

which is the same as in a the case of a smooth feature line associated with a sharpness parameter of 
zero. 

» 

The components of tangent vector weight vector l s can be approximated as l s (s l9 s 2 ) using 
third degree polynomials in sharpness parameters s x and s 2 . The operations performed in generating 
the coefficients for the polynomials correspond to the operations performed in generating the 
coefficients for the polynomials described above for the tangent vector weight vector l c used in 
generating the tangent vector along the smooth feature line (equations (31) through (34)). 
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FIG. 9 depicts a table of coefficients b ti for K=4, 6,. ..,12, with the indexing, in the form of 



the C programming language, b 



K 

T" 2 



[/][/), where "i" is the index of the component of l c , 



taking values in the range / = 0 ? , . . ? 



K 

T +1 



, and "j" is the index of the coefficient by of the 



4 approximating polynomial. 



a 

.u 
05 8 

ao 
S 1 

M 
Q2 

in 



As noted above, the tangent vector weight vectors l c (=l c (s l5 s 2 )) and l s (=l s (s 1 ,s 2 )) (where s, and 
s 2 are the sharpness parameter values at the definition level j=j D for the smooth feature line) are used, 
along with the position d(q) of the vertex v q (0) and the positions d(k) of the neighboring points 
v q (k), k=l,...,K. In particular, the components of the respective tangent vector weight vectors are 
used to weight the positions of the vertex and neighboring points, with the sum comprising the 
tangent vector. That is, the tangent vector e c (q) along the smooth feature line and the tangent vector 
e s (q) across the smooth feature line are generated as 



«c(?) = X \ l c{svS 2 )) c J (i) 

i=0 
i=0 



(37), 



1 3 where ^(s^Sj)^ and l s (s 1 ,s 2 )=l s (where s, and s 2 are the sharpness parameter values at the definition 

14 level j=j D for the smooth feature line), and indices 0,...,K are as shown in FIG. 7. 

15 In some circumstances, it is impractical to use polynomials to generate approximations for 

16 limit point weight vector 1 LP and tangent vector weight vectors l c and l s . For example, if the number 

17 of configurations of smooth feature lines, each requiring a different set of polynomials, exceeds 

1 8 some number, the number of polynomials would increase to an extent that it would be impractical 

19 to use that methodology. An alternative methodology makes use of extrapolation to generate 

20 approximations for limit point weight vector 1 LP and tangent vector weight vectors l c and l s . 
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1 The extrapolation methodology for generating an approximation of limit point weight vector 

2 l L p( s i> s 2) starts the fi rst f ew factors of equation (23) for finite values of index "j." In the 

3 following, it will be assumed that the definition level j D for the smooth feature line is level "zero," 

4 although it will be appreciated that this will not restrict the generality of the description. The matrix 

5 product generated by taking factors j=J,J-l ,...,0 is given by 

6 o 

(38), 

j=J 

N 

Q7 with 

a 

W 8 ^,^(0)^,^) := I K+l (39), 

la 
m 

1^9 where, for the limit case J=0 (in which case there would be no factors in generating the product 

) matrix), the product matrix S scTK ^ LLP (0)(s 1 ,s 2 ) is the identity matrix of size "K+l" by "K+l." A 

r 1 vector of weight values 1 LP for determining the position o(q) of the limit point is obtained as follows: 



W2 



l LP {j)(s x ,s 2 ) = v lp -S scJWP (J)(s x ,s 2 ) (40), 



13 where v LP is as defined above in connection with equations (28) and (29). 

1 4 In one embodiment, a reasonably accurate approximation of limit point weight vector 1 LP can 

15 be obtained using just the vectors corresponding to J=0, 1,2,3 by using polynomial extrapolation. 

16 The convergence of l LP (J)(s 1 ,s 2 ) to l LP (s 1 ,s 2 ) is very rapid, with the error e ( J) = l LP {J) - l LP 

17 approaching zero as 2 if the update rule described above in connection with equations (20) 

18 through (22) are used for the sharpness parameters. The approximation to 1 LP is constructed by 

1 9 generating the third degree polynomial that interpolates the points {x=2~ 2J ,y= 1 LP (J)},5=0, 

20 1,2, 3, and then evaluating this polynomial by extrapolation at the point x=0, that is 
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3 

J=0 



2 where 1 LP (J) is as defined in equation (40) and where the coefficients bj are given by 

-135 



bo 



b, = 



120015 

1270 
120015 

CP < 42 >- 

5 -12192 
S 2 " 120015 

CI _ 131072 

| 3 ~ 120015 

M 

a 

I U 4 This extrapolation methodology can be used for any vertex, whether the vertex is regular or irregular. 

5 In addition, the methodology can be used for any configuration of creases through the respective 

6 vertex, and any value of valence "K" for the respective vertex. 

7 An extrapolation methodology can also be used to generate the tangent vector weight vectors 

8 l c and l s . Generally, in that methodology the matrix products S sc T KL LP (J)(s 1 ,s 2 ) described above 

9 (equations (38) and (39)), which are then 

10 (1) multiplied by a vector v c , if tangent vector weight vector l c is being generated, or a 

1 1 vector v s , if tangent vector weight vector l s is being generated, and 

12 (2) the result from (1) is thereafter multiplied by a dilation factor that maintains the 

1 3 respective tangent vector weight vector in normalized form. 

14 The approximations of the tangent vector weight vectors l c and l s are generated from 
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l c (j)(s lf s 2 ) = d{K) J v c S scTKLLP (j)(s, y s 2 ) 
l s {j){s x ,s 2 ) = d(K) J 



(43), 



where d(K) represents the dilation factor. The vectors v c and v s are given by 



2*(0) 2x{\) 
0, cos — — , cos — — , . . . , cos 



2x(K- 1) 



K 



K 



K ) 



0, sin 



2*(0) . 2*(l) 



, sin- 



K 



sin - 



(44), 



and the dilation factor d(K) is the reciprocal of one of the eigenvalues of the subdivision matrix for 
the smooth case, namely, 



d(K) = 



1 



3 1 In 
— + — cos — 
8 4 K 



(45). 



Approximations to the tangent vector weight vectors l c and l s are given by 



J=0 



(46), 



where the coefficients bj have the form as given in equation (42). This extrapolation methodology 
can be used in connection with all configurations. The tangent vector generated using the weight 
vector l c will be oriented approximately in the direction of the edge [1,0] between vertex v q (0) and 
point v q (l), and the tangent vector generated using the weight vector l s will be oriented 
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1 approximately perpendicular to the edge [1,0]. Although the angle between edges [ 1 ,0] and [0,L+ 1 ] 

2 for a given level "j" in a subdivision surface is typically not one hundred and eighty degrees, as the 

3 level "j " approaches infinity, the angle will approach one hundred and eighty degrees. Accordingly, 

4 in the limit, that is, as "j" approaches infinity, the tangent vector along the edge on one side of the 

5 limit point associated with a vertex v q (0) will approach the opposite direction as the tangent vector 

6 along the edge on the other side of the limit point associated with the same vertex. Similarly, the 

7 tangent vector across the edge on one side of the limit point associated with the vertex v q (0) will 

8 approach the opposite direction as the tangent vector across the edge on the other side of the limit 
p 9 point associated with the same vertex. 

a 

IJIO With this background, operations performed by the computer graphics system 10 in 

PJ 

I j 1 generating the weight vectors 1 LP , l c and l s useful in generating the limit point and tangent vectors 

*£l2 for a subdivision surface for a triangular mesh will be described in connection with the flowchart 

ru 

n 13 depicted in FIG. 10. With reference to that FIG., the computer graphics system 10 first initializes 

4 the subdivision matrices S sc T ICL (s 1 (j),s 2 (j)), j=0,l ,2 (step 100) as described above in connection with 

Gl5 equation (19). Thereafter, the computer graphics system 10 generates the matrix products 

M 

Ql6 S SCJ ^ LP (J)(s l3 s 2 ) for J-2,3 (step 101). For J=l, S scXKjULP (l)(s 1 ,s 2 )=S SC)TXL (s 1 ,s 2 ), and, as indicated 

n ^l 7 above, for J=0, S SCJXL>LP (0)(s 1 ,s 2 ) is the "K+r-by-"K+r identity matrix. 

1 8 After the computer graphics system 10 has generated the matrix products S sc T K L LP , it takes 

19 the first row of each matrix product S sc T I ^ L LP (J)(s 1 ,s 2 ), J=0, 1, 2, 3, and generates component-wise 

20 an approximation to limit point weight vector 1 LP in accordance with the extrapolation formula in 

21 equation (41) (step 102). In addition, the computer graphics system 10 uses the matrix products 

22 S scT K L LP (J)(s 1 ,s 2) , dilation factor d(K) and vectors v c and v s to generate the vectors 1 C (J) and 1 S (J), 

23 J=l , 2, and 3 as described above in connection with equation (43) (step 1 03). For J=0, the respective 

24 weight vectors are l c (0)=v c and l s (0)=v s . Following step 103, the computer graphics system 10 can 

25 use 1 C (J) and 1 S (J) to generate the approximations to the tangent vector weight vectors l c and l s in 

26 accordance with equation (46) (step 1 04). After the limit point weight vector 1 LP and tangent vector 

27 weight vectors l c and l s have been generated, the limit point and tangent vectors can be generated 

28 as described above (reference equations (30) and (37), respectively), and the normal vector can also 

29 be generated as the cross product between the tangent vectors (step 105). 
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1 Operations performed by the computer graphics system 10 in connection with generating a 

2 smooth feature line in the case of a subdivision surface defined by a quadrilateral mesh will be 

3 described in connection with FIGS. 1 1 through 13 and the flow chart in FIG. 10. With reference 

4 to FIG. 11, that FIG. depicts a portion of a subdivision surface defined by a triangular mesh with 

5 vertex v q (vertex v q =v q (0)) and points v q ( 1 ) through v q (2K) with indices numbered as depicted in the 

6 FIG. It will be appreciated that points v q ( 1 ) through v q (K) are first-order neighbors of vertex v q , and 

7 points v q (K+ 1 ) through v q (2K) are second-order neighbors of vertex v q . The vertices are labeled such 

8 that a crease passes through point v q ( 1 ), vertex v q (0) and point v q (L+ 1 ) along edges v q ( 1 ) ,v q (which 

La 

L9 will be referred to as edge (1,0)) and v q ,v q (L+l) (which will be referred to as edge (0,L+1)). The 

□ 0 value of the sharpness parameter associated with the edge ( 1 ,0) will be referred to as s j and the value 

01 

1 of the sharpness parameter associated with the edge (0,L+1) will be referred to as s 2 . Generally the 

^ 2 subdivision rule described in equations (13) and (14) can be used to determine the positions of 

FU3 respective vertices, using, instead of the subdivision matrix described above in connection with 

s 

J 4 equation (15), the subdivision matrix 
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8 Sj 



Ms*) 

3/ x 1 

3 
8 



J_ 
16 



4 



s 3 if I = 0, m = 0 

z/ / = 0, w = 1 or L + 1 

i//= 0,m = 2,...,Z, 
or I = 0, m = L + 2,...,K 

ifl=0,m=K + l,...,2K 

if I = l, m = 0 or 1 
if I = L + l,m= 2, K, 
K+lor2K 

ifl=L+l,m=0orL+l 
if 1= L+l,m= L,L+2, 
K+ lor K+ L+ 1 

i//= 2,...,L,m = 0 
or /= 1+ 2,...,K,m = 0 
or / = w = 2,...,Z, 
or I = m = L+ 2,...K 

ifl= 2,...,L,m= I- 1 
or / = 2,...,L, m = I + 1 
or /= 2,...,L,m= K+ I- 1 
or / = 2,...,L,m= K+ I 
or 1= L+ 2,...,K,m= I- 1 
or / = L + 2,...,K- 1, /w = /+ 1 
or / = K,m- 1 

or/= L+2,...,K,m= K+l- 1 

or /= L+ 2,...,K,m= K+ I 
if 1= K+ l,...£K- l,m= 0, 

l-K,l- K+lorl 
orl=2K,m= 0, K, 1 or 2/: 

otherwise 



(47), 
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is used, where, as with equation (19), s 3 = ~\S\ + s 2 j ,the average of the values of the parameters 



s, and s 2 associated with the edges (1,0) and (0,L+1). 

As with the case of a subdivision surface defined by a triangular mesh, to ensure that the 
direction of the respective tangent vectors vary continuously at and near irregular vertices at the 
various subdivision levels, the value(s) of the sharpness parameter(s) are gradually reduced as the 
subdivision level of the mesh increases. In one embodiment, the computer graphics system varies 
the value(s) of the sharpness parameter(s) in the case of a subdivision surface defined by a 
quadrilateral mesh in the same manner as in the case of a subdivision surface defined by a triangular 
mesh (reference equations (20) through (22) above). 

As with the case of a subdivision surface defined by a triangular mesh, to determine the 
positions of the limit point associated with a vertex v q (0) in the case of a subdivision surface defined 
by a quadrilateral mesh, the infinite subdivision matrix product 



is evaluated, where S scQKL (s 1 (j) 5 s 2 (j)) is the subdivision matrix from equation (47) for sharpness 
parameter values s, and s 2 associated with the j-th level and index j D refers to the definition level of 
the smooth crease. In the left-hand side of equation (48), for the arguments s 1 and s 2 , s^s,^) and 
s 2 =s 2 (j D ). As in the case of a subdivision surface defined by a triangular mesh, the matrix product 
in equation (48) converges to a matrix with identical rows, in this case each row being of size 2K+1 . 
A vector of weight values 1 LP for determining the position a(q) of a limit point associated with vertex 
v q (0) is obtained as follows: 



where vector v LP is a vector of weights that are used in determining the position of a limit point of 
a subdivision surface defined by a quadrilateral mesh in the absence of a smooth feature line, that 
is, 





(49), 



v 



LP 




(50), 
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1 where, in equation (50), there are "K" components having the value "four" and "K" components with 

2 value "one." 

3 The components l LP (s 1 ,s 2 ) i , i=0,...,2K, of the row vector l LP (s„s 2 ), and the positions d(0) of 

4 the vertex v q (0) and d(i) (i=l ,...,2K) of the neighboring points v q (i) are used to determine the position 

5 a(q) of a limit point associated with the vertex v q (0) as follows 

*M = EM^Jp) (51) ' 

S Ql where, as suggested above, i equ N e (cj[,j)u N f {q 9 j) (reference FIG. 11). 

~ 1 8 As with the case of a subdivision surface defined by a triangular mesh, the components of 

?y - 

M 9 the limit point weight vector l LP (s,,s 2 ) can be approximated using polynomial expansions. The 

*J5 

pJO components of limit point weight vector l LP (s l3 s 2 ) can be satisfactorily approximated by second 

y 1 degree polynomials in s, and s 2 using six values per component and per valence. If the value of "K" 

ru K 

UL2 is even, and if L - — , symmetry considerations allow the number of values per component to be 

M ^ 

CD 

1 \ 3 reduced. In that case, the components of limit point weight vector 1 LP satisfy the following symmetry 

14 relations: 

15 1. the "i-th" component of limit point weight vector l L p(s„S2) equals the "K+2-i-th" 

16 component of limit point weight vector l LP (s 1 ,s 2 ), that is, 

1 7 (l Lp (s x , s 2 )) . = (l LP (s x , s 2 )) K+2 i ,for/=2,...,y,andthe "i-th" component of 

18 limit point weight vector l LP (Sj,s 2 ) equals the "3K+l-i-th" component of limit point 

19 weight vector l LP (s„s 2 ), that is, (hp{ s \ 9 s 2)). = {hp{ s \^ s i)) 3K+l _^ for 

K 

20 /= K+ l,...,K+— ; 
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2. (l Lp [s x , s 2 )) = (/i/>(5,,5 2 ))^ +2 for 1=1,...,—+!, and 



(l LP {s l ,s 2 )) j = (/L/»(-s,,5 2 )) 2 ^ +£+i _.,for i = K+ l,...,K + y + 1; 



3. (/ £ p(^,5 2 )) 0 = (/^(5 2 ,5,)) o ; and 

4. (^(-Si^ 2 ))£ +1 = (/ z/ ,(5 2 ,5,)) £+i , if "K" is a multiple of "four," and 

4 4 

M*i>* 2 )L^j +1 = M-^-OL^, ' if " K " is not a muitipie ° f ,,four - M 



Given the symmetry relations described immediately above, it is only necessary to generate 

K 



approximations for limit point weight vector 1 LP for j=0,. 



+ 1 and 



4 



+ 1 , and, in addition, symmetric polynomials can be used for i = 0 for any 



even value of "K," as well as for / = — + 1 if "K" is a multiple of "four," or / = K + 



4 



+ 1 if 



"K" is even but not a multiple of "four." 

As in the case of a subdivision surface defined by a triangular mesh, in one embodiment i the 

K 

case of a subdivision surface defined by a quadrilateral mesh, regardless of whether L = — , the 



approximating polynomials are computed using a least-squares Chebyshev approximation 
methodology (reference equation (27) above). In this case, the coefficients by (j=0,...,3 in the 
symmetric case, or j=0,...,5 in the asymmetric case) for the polynomial 
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(>4 z 6 .o+^i+^) + M s i 2+5 2 2 ) + te ( 52 >> 

in the symmetric case, or the polynomial 

(lu), * *>o + ^fl + b i2 s 2 + b n s 2 x + b iA s x s 2 + ft /5 5 2 2 (53), 

in the asymmetric case, are determined by a least squares methodology. FIG. 12 depicts a table of 
coefficients b^ for K=4, L=2, for values of "i" i=0; 1; 2; 5. After the coefficients b^ have been 
generated (equations (52) and (53)), they are used to generate the limit point weight vector 1 LP in the 
same manner as in the case of a subdivision surface defined by a triangular mesh, which can be used 
to generate the position of the limit point a(q) as described above in connection with equation (51). 

The tangent vectors in the case of a subdivision surface defined by a quadrilateral mesh are 
generated in a manner similar to that used to generate tangent vectors in the case of a subdivision 
surface defined by a triangular mesh. Generally, in that methodology the matrix products 
S sc Q K L LP (J)(s 1 ,s 2 ) as described above (equation (48)) are generated, which are then 

(1) multiplied by a vector v c , if tangent vector weight vector l c is being generated, or a 
vector v s , if tangent vector weight vector l s is being generated, and 

(2) the result from (1) is thereafter multiplied by a dilation factor that maintains the 
respective tangent vector in normalized form, 

as in the case of the tangent vector weight vectors in the case of a triangular mesh. Initially, the 
vectors l c (J)(s 1? s 2 ) and l s (J)(s 1? s 2 ) are generated as 

l c (s l9 s 2 ) = d{K) J v c ^S SC9QiKLtLP (j)(s l9 s 2 ) 

(54). 

l s {s l9 s 2 ) = d(K) J 

In the case of a subdivision surface defined by a quadrilateral mesh, the respective vectors v c and 
v s are given by 
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131 

ru 

m 

8 i 

o 

fy 2 
Q 

M3 



(vcl = 



M, - 



0 



COS 



2x(i-l) 



ifi=0 
if i = 1,..., AT 



2*(/-iE-l) . 2x(i-K) 
cos — +cos — if i = K + l,...,2K 



K 



K 



0 



A K sin 



2n{i-\) 
K 



ifi=0 

if i = 1,. . . , K 



2x(i-K-\) . 2x(i-K) 
sin — + sin — if i = K + 1,...,2a 



(55), 



and the dilation factor is given as the reciprocal of the subdominant eigenvalue of the subdivision 
matrix (equation 47) for a smooth interior vertex: 

1 

A K 1 (56), 
— + - 
16 4 



d(K) = 



4 

5 
6 



where A K is defined by 



= 1 + CO; 



(2k\ 
a — 
V K) 



+ co{^) 



2 9+ cos- 



2^ 
A 7 ) 



(57). 



7 The tangent vector weight vectors l c and l s are used to generate the tangent vectors as 

8 follows: 
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(58), 



where the indexing is as shown in FIG. 11. 

As in the case of a subdivision surface defined by a triangular mesh, in the case of a 
subdivision surface defined by a quadrilateral mesh the tangent vector e c (q) generated for a vertex 
v q (0) generated using the tangent vector weight vector l c will be oriented approximately along the 
smooth feature line through the limit point associated with the vertex v q (0), and the tangent vector 
e s (q) that is generated using the tangent vector weight vector l s will be oriented approximately 
perpendicular to the smooth feature line through the limit point associated with the vertex v q (0). 

Also, as with the case of a subdivision surface defined by a triangular mesh, the tangent 
vector weight vectors in the case of a subdivision surface defined by a quadrilateral mesh can be 
approximated using polynomial approximations. The components of tangent vector weight vectors 
l c (s ls s 2 ) and l s (s 1 ,s 2 ) can be satisfactorily approximated by third degree polynomials in s, and s 2 
using, in general, ten values per component and per valence If the value of "K" is even, and 

K 

ifL= — , symmetry considerations allow the number of values per component to be reduced. In 

that case, the components of tangent vector weight vector l c satisfy the following symmetry 
relations: 

1 . the "i-th" component of tangent vector weight vector l c (s l5 s 2 ) equals the "K+2-i-th" 
component of tangent vector weight vector tangent vector weight vector l c (s, ,s 2 ), that 




, for i = 2, . . . , — , and the "i-th" component 



of tangent vector weight vector l c (s l5 s 2 ) equals the "3K+1 -i-th" component of tangent 
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vector weight vector lc(s„s 2 ), that is, (/ c (j,,5 2 )) = (/ c (s,,S 2 )) 3 ^ +i . , for 



/ = a:+ 1,...,^"+ — ; 



2. (/cfci^i)),. = " ('cfci'^))*^.,. f 0 r * = 1, and 



(/ c (j„5 2 )) = " ('cfo^L * . ,.> for * = ^+ 1,...,^+ y + 1; 

3. (/ c (5,,5 2 )) o = - (/ c (j 2 ,5,)) o ; and 

4. (/c( s i'' s 2))^ = ~ (^c( 5 2» ,s i))£ ' if " K " is a multiple 0I " "four," and 
(/ c (5 1 ,5 2 ))^ + ^j +i = - (/ c (-S 2 ,s,))^ + |Kj +i » if " K " is not a multiple of "four." 



Given the symmetry relations described immediately above, it is only necessary to generate 



approximations for i - 0,. . . . 



K 



+ 1 and /= K+ \,...,K + 



K 



+ 1 , and, in addition, anti- 



K 

symmetric polynomials can be used for i = 0 for any even value of "K," as well as for / = — +1 

4 



if "K" is a multiple of "four," or i - K + 



K_ 
4 



+ 1 if "K" is even but not a multiple of "four." 



As in the case of a subdivision surface defined by a triangular mesh, in one embodiment, the 
approximating polynomials are also computed using a least squares Chebyshev approximation 
methodology similar to that used in connection with the limit point approximation described above. 
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1 The matrix product of equation (48) is evaluated at points of a grid in the (s 1? s 2 ) region made up of 

2 roots of Chebyshev polynomials (equation (27)). In this case, the coefficients (j=0,...,3, in the 

3 anti-symmetric case, or j=0,...,9, in the non-symmetric case) for the polynomial 

4 (U • ^(st-sJ+bnfc-^+bnfa (59), 

5 in the anti-symmetric case, or the polynomial 

(oil), 

u b iS sl + b i6 sl + b n s 2 x s 2 + b i% s x 4 + b i9 s\ 

s 

£57 in the non-symmetric case, are determined using the least squares methodology. Operations 

01 

py 8 performed in connection with the least squares methodology correspond to those described above 

s ~! 9 in the case of a subdivision surface defined by a triangular mesh. The components of tangent vector 

WO weight vector l s can be approximated as l s (s,,s 2 ) using third degree polynomials in sharpness 

j^L 1 parameters s, and s 2 . The operations performed in generating the coefficients for the polynomials 

Hjl 2 correspond to the operations performed in generating the coefficients for the polynomials described 

y\ 3 above for the tangent vector weight vector l c used in generating the tangent vector along the smooth 

Q 

~14 feature line. FIG. 1 3 depicts a table of coefficients b^ for K=4, L=2, corresponding to i=0; 1 ; 2; 5. 

15 After the tangent vector weight vectors l c and l s have been generated, they are used to 

16 generate the tangent vectors e c and e s as described above in connection with equation (58). 

17 As in the case of a subdivision surface defined by a triangular mesh, in some circumstances, 

18 it is impractical to use polynomials to generate approximations for limit point weight vector 1 LP and 

19 tangent vector weight vectors l c and l s . An alternative methodology makes use of polynomial 

20 extrapolation to generate approximations for limit point weight vector 1 LP and tangent vector weight 

2 1 vectors l c and l s . This extrapolation methodology is similar to the methodology described above in 

22 connection with equations 38-46. 

23 Operations performed by the computer graphics system 1 0 in connection with generating the 

24 limit point weight vector 1 LP and tangent vector weight vectors l c and l s for a subdivision surface 

25 defined by a quadrilateral mesh, using the polynomial extrapolation methodology, are similar to 

26 those described above in connection with FIG. 10 for a subdivision surface defined by a triangular 

27 mesh, and reference should be made to FIG. 10 for that description. 
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1 The arrangement provides a number of advantages. In particular, the invention provides an 

2 arrangement for generating a subdivision surface from a representation defined by a control mesh 

3 at a selected level, by applying the same parameterized subdivision rule at each level, with the 

4 value(s) of the parameter(s) varying at at least some of the levels. The invention specifically 

5 provides for the generation of a smooth feature line in a subdivision surface using the same 

6 parameterized subdivision rule at each level, with the value(s) of the parameter(s) defining the 

7 smooth feature line varying at at least some of the levels. The resulting subdivision surface is of 

8 relatively high quality and can be efficiently evaluated. 

g9 It will be appreciated that numerous changes and modifications may be made to the 

Q 0 arrangement as described herein. For example, it will be appreciated that, although the arrangement 

1 has been described in connection with meshes of triangular and quadrangular faces, it will be 

jjj 2 appreciated that the arrangement may find utility in connection with meshes having faces of different 

fU3 polygonal structures. 

F*jL4 Furthermore, although specific subdivision rules (equations (19) and (46)) have been 

Hjjl 5 described, it will be appreciated that other subdivision rules may find utility. 

Ml 6 In addition, although a methodology has been described for updating the sharpness 

a 

~jl7 parameters s, and s 2 , and although various methodologies have been described for approximating 

18 weight vectors 1 LP , l c and l s , it will be appreciated that other parameter update and approximation 

1 9 methodologies may find utility. 

20 Furthermore, in connection with the extrapolation methodologies descrbied herein, it will be 

21 appreciated that the number of coefficients b j5 and the upper limit ("three") on the summations in 

22 equations (41) and (46) may differ if a polynomial other than of order "three" is used in the 

23 extrapolation. 

24 It will be appreciated that a system in accordance with the invention can be constructed in 

25 whole or in part from special purpose hardware or a general purpose computer system, or any 

26 combination thereof, any portion of which may be controlled by a suitable program. Any program 

27 may in whole or in part comprise part of or be stored on the system in a conventional manner, or it 

28 may in whole or in part be provided in to the system over a network or other mechanism for 

29 transferring information in a conventional manner. In addition, it will be appreciated that the system 

30 may be operated and/or otherwise controlled by means of information provided by an operator using 
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operator input elements (not shown) which may be connected directly to the system or which may 
transfer the information to the system over a network or other mechanism for transferring 
information in a conventional manner. 

The foregoing description has been limited to a specific embodiment of this invention. It will 
be apparent, however, that various variations and modifications may be made to the invention, with 
the attainment of some or all of the advantages of the invention. It is the object of the appended 
claims to cover these and such other variations and modifications as come within the true spirit and 
scope of the invention. 

What is claimed as new and desired to be secured by Letters Patent of the United States is: 
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